Exchange coupling in semiconductor nanostructures: 
Validity and limitations of the Heitler-London approach 
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The exchange coupling of the spins of two electrons in double well potentials in a semiconductor 
background is calculated within the Heitler-London (HL) approximation. Atomic and quantum dot 
types of confining potentials are considered, and a systematic analysis for the source of inaccuracies 
in the HL approach is presented. For the strongly confining coulombic atomic potentials in the H2 
molecule, the most dramatic failure occurs at very large interatomic distances, where HL predicts a 
triplet ground state, both in 3D and in 2D, coming from the absence of electron-electron correlation 
effects in this approach. For a 2D double well potential, failures are identified at relatively smaller 
interdot distances, and may be attributed to the less confining nature of the potential, leading 
to larger overlap. We find that in the double dot case, the range of validity of HL is improved 
(restricted) in a related 3D (ID) model, and that results always tend to become more reliable as the 
interdot distance increases. Our analysis of the exchange coupling is of relevance to the exchange 
gate quantum computer architectures in semiconductors. 

PACS numbers: 03.67.Lx, 85.35. Be, 73.21. La, 85.30.-z 



I. INTRODUCTION 

Semiconductor spin-based quantum computation 
emerged from theoretical proposals showing that the 
required universal gatesi could be implemented through 
physical operations involving single electrons bound to 
an array of quantum dots2^i^*SiLSi2iiS or to donorsiii 
The localized spin of each electron serves as the single 
qubit by virtue of the 2-level spin-dynamics, and the 
(electrostatic) quantum mechanical exchange coupling 
between the two electrons can be used to entangle the 
qubits through the "exchange gate" . 

Accurate calculations of exchange coupling in 
molecules (or artificial molecules) are extremely dif- 
ficult requiring numerically intensive self-consistent 
solutions of Hartree-Fock equations. In semiconductor 
nanostructures, particularly in the context of quan- 
tum computer architectures, such intensive quantum 
chemistry-type numerical calculations are unwarranted 
for several reasons. First, the parameters entering such 
calculations are essentially only very approximately 
known in semiconductor nanostructures and, therefore, 
extremely accurate quantitative calculations for the 
exchange energy in these nanostructures is not partic- 
ularly meaningful. Second, the quantum chemical type 
exchange calculations are not particularly valid in the 
semiconductor environment. 

One of the simplest and most successful methods for 
the calculation of exchange coupling in small molecules 
was proposed in 1927 by Heitler and LondoniiSsi 3 . The 
basic assumption of the Heitler-London (HL) method is 
that the many-electron lowest energy wavefunctions in 
molecules may be written in terms of the one-electron 
ground-state orbitals of the isolated constituent atoms. 



This immediately requires that the overlap (S) among 
neighboring orbitals be small, a condition not always ful- 
filled in real molecules. In fact, for the equilibrium in- 
teratomic distance in the H2 molecule (1.5 as, with as 
representing the Bohr radius), the overlap is larger than 
0.7, and the HL ground-state energy is overestimated by 
~ 1.5 eV, while for interatomic distances R > 5 0,3 the 
overlap becomes smaller than 0.1 and an excellent agree- 
ment is obtained between HL and the experimentally ob- 
served ground-state energy^ The HL approximation is 
thus expected to be valid for S <C 1, or equivalently, for 
well separated atoms. 

A somewhat surprising anomaly of the HL method ap- 
plied to H2 is that it predicts a triplet ground-state for 
R > 50 as, in contradiction with the well established 
result that, for a spin-independent Hamiltonian and in 
the absence of magnetic fields, the ground-state of an 
even-electron system in a symmetric potential must be a 
singlet. 15 This anomaly in HL comes from the singular- 
nature of the Coulomb potential, 13 leading to a logarith- 
mic dependence of the electron-electron repulsion energy 
with R which becomes the dominant term at large R, 
resulting in a negative value of the exchange coupling 
J(R) = Etnpiet-E S i n gi e t in this limit. An asymptotically 
correct expression for J(R) for the hydrogen molecule 
was obtained by Herring and Flicker^ who have also 
shown that the HL results for J(R) are in excellent quan- 
titative agreement with the asymptotic expression up to 
R values very near the crossover to the unphysical nega- 
tive range of J(R). 

Many calculations of exchange coupling in semicon- 
ductor quantum dots are based in the HL approximation 
applied to 2D model potentials. It has been pointed out 
that, for some parameter values of a quartic model poten- 
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tial proposed for the study of 2D double quantum dots, 
unphysical (negative) values of J are obtained within the 
HL approximation^ This failure of HL may be due to 
the reduced dimensionality of the potential, since the 
electron-electron Coulomb interaction in 2D is expected 
to be larger than in 3D, or to some intrinsic limitation of 
this approach when applied to weakly confining gated po- 
tentials as compared to atomic- like molecular potentials. 
We address these questions here by considering the 2D 
hydrogen molecule and the two-electron double quantum 
dot quartic potential within HL, including its modified 
versions in 3D and ID. 

Our study is complementary to detailed and numer- 
ically intensive calculations of exchange. We find that 
the reliability of HL depends primarily on the form of the 
confining potential: (i) For weakly confining (low barrier) 
potentials, the range of validity is generally improved in 
larger space dimensions; (ii) For strongly confining po- 
tentials, illustrated here by the coulombic potential, the 
space dimensionality is not a decisive parameter; (iii) Un- 
reliable and even unphysical results may be attributed 
to non-additive model potentials intended to simulate a 
double- well environment. 

This paper is organized as follows: In Sec.[H]we quickly 
review the HL method for the H 2 molecule, and consider 
the two-dimensional analogue of this problem; in Sec. lIIII 
we present the HL solution for the two-electron problem 
in the 2D double quantum dot quartic potential intro- 
duced in Ref. and compare it to its modified versions 
in 3D and in ID. Further discussions and conclusions are 
presented in Sec. IIVI 



II. THE HL METHOD FOR THE H 2 
MOLECULE IN 2D AND 3D 



The Hamiltonian of the hydrogen molecule with nu- 
clear coordinates Ra and R# is written 



H = T 1 + T 2 + V(r 1 ) + V(r 2 ) + 



e 

ri2 



e 
~R 



(1) 



where V(r 4 ) = Va^) + V B (r t ) = -f- - f~ b , r ia = \n - 
R-a|, — |r$ — Rb|, Yi (i = 1, 2) are the electronic 
coordinates, Ti is the kinetic energy operator for electron 
i, ri2 = |ri — r2 1 and R — |Ra — Rb| is the inter- nuclear 
separation. Starting from normalized hydrogen atomic 
orbitals centered at Ra and Rg, which we denote for 
each electron i as a(i) and b(i), the HL lowest singlet 
(+) and triplet (— ) 2-electron states (we omit the spin- 
dependent part here) are written as 



l±) = 



1 



x/2(l±S 2 ) 



[o(l)ft(2)±6(l)o(2)], (2) 



with S = (a(i)\b(i)) giving the overlap integral. It is 
convenient to cast the exchange coupling, J — E tr ipiet — 



Esingiet = (-|-ffj-) - (+\H\+) in the following form 
2S 2 



J = 



1-S 4 



(W-C) 



(3) 



where 



W = (a(l)b(2)\v\a(l)b(2) - —a(2)b(l)) 
with v = V(ri) + V(r 2 ) - V B (r 2 ) - V A (ri), and 



(4) 



C=(a(l)6(2)| — |-o(l)6(2) 
ri2 



S 2 



a(2)6(l)> . (5) 



The two terms W and C are always positive, and corre- 
spond to energy terms which contribute differently to the 
singlet and triplet states: W is related to the "covalent" 
energy, which favors the singlet state providing molecu- 
lar bonding, while C is related to the electron-electron 
Coulomb repulsion, favoring the triplet state. Note that 
both terms are independent of the kinetic energy. As 
discussed in the introduction, physically J > 0, therefore 
we expect W > C. 

The HL solution was presented long ago for the 3D 
caseji^ where an isolated atomic orbital has the form 
a3D(i) = e~ Tia /y/ir , and equivalently for 630(1)- Dis- 
tances are given in units of a B = h 2 /(me 2 ). The 3D 
overlap is S 3D = exp(-ii) (l + R + R 2 /3) . 

The 2D version of this problem has not been previ- 
ously discussed within HL. We present it here since it is 
useful in addressing the effects of dimensionality, particu- 
larly regarding the electron-electron Coulomb potential. 
The ground-state orbital for the 2D hydrogenic atom 
i o 17 ' 18 a 2 o(i) = 4exp(— 2r ia )/y/2n , leading to the over- 
lap S 2D = 2i? 2 [A' (2i?) + A' 1 (2i?)/i?] , where K and 
Ki are modified Bessel functions of the second kind. The 
problem in 2D is solved following the same procedure as 
in 3D, in particular using a 2D version of the spheroidal 
coordinates A = (r a + rh) / R and [i = (r a — rb)/R de- 
scribed in detail in Ref. [lj. Contrary to the 3D case, 
for which analytic expressions are obtained for all terms, 
in 2D the electron-electron Coulomb interaction terms 
contributing to C in Eq. JSJ are obtained numerically. 

In Fig. ^ w e present the results for the overlap (£>), the 
covalent energy (W), the Coulomb energy (C), and the 
exchange (J) in two and three dimensions. Energies are 
given in units of Ry= me 4 /(2h 2 ). The first thing that 
shows up is that, although the overlap S 2 d(R) is smaller 
than S3d(R) for all R, both W and C contributions, as 
well as their difference W — C, are larger in 2D than in 
3D, which reflects the enhanced Coulomb interaction due 
to the electronic charge additional confinement in 2D as 
compared to 3D. Fig.^c) illustrates the well known arti- 
fact of HL for the 3D hydrogen problemi& giving negative 
J at very large distances (above R ~ 50as). At shorter 
distances, the comparison of J given by Eq. JHJl with the 
most accurate asymptotic results available in the liter- 
ature for 2D and SDiSiii is very good for a wide range 
of distances. Within our numerical accuracy, in 2D the 
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FIG. 1: (Color online) (a) Overlap between the hydrogen 
ground-state orbitals centered at nuclear sites a distance R 
apart in 2D and 3D. (b) Covalent (W) and Coulomb (C) en- 
ergy terms contributing to the electronic exchange within the 
HL approach in the hydrogen molecule in 2D and 3D, in units 
of Ry. The downturn region at short distances (R < 4) is less 
reliable due to the large overlap, (c) The solid lines are the 
asymptotic values for the electronic exchange for the hydrogen 
molecule in 2D and 3D^ii j^ymptotic = 15 2R v/4 exp( _ 47? ) 

and J^ m P totic = 0.8i? 5/2 exp(-27?), that we compare here 
with the Heitler-London results (dashed lines) . The compari- 
son is very good up to the large distances where HL fails. The 
HL results are dominated by the prefactor S 2 , see Eq. (J^J , ex- 
cept for the negative- J anomaly, indicated here by the steep 
drop in the HL curves. 



critical R at which J becomes negative is between 40 and 
lOOas, comparable to the 3D value. At very short dis- 
tances, HL becomes less reliable due to the large overlap 
[see Fig. [T][a)]. Arbitrarily denning a cutoff at S = 0.1, 
we estimate the minimum distance at which we can trust 
the HL results to be R^ n ~ 5as, and R^ n ~ 2.5as, 
since the overlap for the 2D system decays much faster 
with R than for 3D. 



III. THE HL METHOD FOR THE QUARTIC 
DOUBLE-DOT POTENTIAL 

Quantum dots are defined by gate-generated poten- 
tials modifying a two-dimensional electron gas (2DEG) 
environment in semiconductor heterostructures. Each 
quantum dot is usually modeled by a harmonic well with 

V 3 (r) = ^i(r - Kj) 2 where Rj = +R/2x, -R/2x 
for j = A.B respectively. Following Burkard et al.fi the 
coupling of the dots in a 2D system is modeled by the 
quartic potential 



V(x,y) 



TflLOZ, 



(If 



R 2 



(6) 



which satisfies V(x w ±i?/2,y) = Vj(r). Note that the 
barrier height increases with R. The Hamiltonian for two 
electrons in a double quantum dot can be written as 



H = T 1 +T 2 + V(r 1 ) + V(r 2 ) 



£7-12 



(7) 



where e is the semiconductor dielectric constant. For the 
results presented below, distances are given in units of 
a* = h 2 e/(m*e 2 ) and energies in Ry* = m*e i /(2h 2 e 2 ), 
where m* is the effective mass in the host semiconductor. 
Note that, contrary to the H 2 case, V ^ Va + Vb- The 
single electron Hamiltonian is 

h i = T i +V A (r i ) , (8) 

the corresponding ground state orbital is a 2 D(i) = 
Jj- exp [(a* ± f ) 2 + 2/f] }, with (3 = y/m*u /h, 

and the overlap is given by S 2 d = exp(— /3 2 (i?/2) 2 ). 

We point out that the physical significance of the co- 
valent term W defined in Eq. Q is somewhat differ- 
ent here as compared to the H2 system. From Eq. Q, 
the operator involved in W is an effective potential v = 
V(ri) + V(r 2 ) — Vs(r 2 ) — Vk(ri) which, for H 2 , results in 
v = Vb(i*i) + Va{y 2 ), related to the energy of an "A-atom 
electron" due to the potential Vb of the £?-atom and vice 
versa. In the double quantum dot, v ^ Vs(ri) + Va(t 2 ), 
instead the effective potential v involves differences be- 
tween the quartic potential and quadratic terms. In this 
way, the "electron in the A-dot" does not "feel" the po- 
tential Vb, but an approximated version of it. 

Using the expressions for W and C given in Eqs. Q 
and © we get fi 



3, / /3 2 R 2 
W 2D = -hajQ I 1 + — ^— 



and 



C 2 D = ftWQ C 



1 



(3 2 R 2 



(9) 



(10) 



where Iq is the modified Bessel function of the first kind 
and 



7T e 2 p/e 
2 Tiluq 



(11) 



is a parameter introduced in Ref. |3J as the ratio between 
the Coulomb and the confining energy. Note that c is 
also completely defined by ujq, since (3 ~ a/cjq. 

In order to investigate the effect of dimension- 
ality, it is instructive to compare the 2D system 
with a similar 3D problem. From the 2D po- 
tential in Eq. ©, we define a 3D counterpart as 

2 

Vzd{x, y, z) = V(x, y) + ^p-z 2 , with the ground state or- 
bital a 3D {i) = f^exp|-/3 2 (xi ± f ) 2 
leading to S 3D — S 2D , W 3D = W 2D and 



Vi 



/»}. 



Q 



3D 



-hujQ c 



1 



IT Erf(/3i?/\/2) 

2 ~jm 



(12) 
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5 10 15 20 25 30 35 
R(a*) 

FIG. 2: (Color online) (a) Overlap (S = Sid = S 2D = 
S3d) for three different values of the parameter c [defined 
in Eq. (lilt ] for the double quantum dot problem, (b), (c), 
and (d) Covalent (W = Wid = W2D = Ww) and Coulomb 
(C) energy terms for the double quantum dot problem in ID, 
2D, and 3D. When the overlap is large, HL results are not 
reliable and, for larger values of c, lead to negative exchange 
J. 



where Erf (x) — -4= J* e u ' 2 du is the Error Function. 



For completeness, we also discuss the ID problem. 

2 



The model potential is V\d (x) : 
The single-electron ground 
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X 



(f) 

a w (i 



/R 2 



state is 

■exp|-/3 2 (xi± §) 2 /2j, leading to Sid = S 2 d, 

W\d — W 2 D- The Coulomb term corresponds to the 
e = value of the function 

e -0 2 x 2 /2 



CiD(e) = -hcj c 

7T 



dx ■ 



x 



- Uujo c S 2 



dx ■ 



-0 2 x 2 /2 cosh ( / 32 Ea .) 



(13) 



where x = \x\ — xa|- I n l ne e — * limit, the singular 
nature of the Coulomb potential in ID systems leads to 
a logarithmic divergence of this expression. We overcome 
this problem phenomenologically by introducing a cutoff 
distance (the results presented in Fig. correspond to 
e = 0.2a*) to simulate the correlation effects that would 
lead to a lower probability density at x\ = x%- Such 
correlation effects are not included in HLpi^ 

Fig. 13 shows the behavior of S, W and C vs R for 
different values of c and dimensionalities. The HL ap- 
proach breaks down for small R for a certain value of 



c which depends on dimensionality. For 2D, the break- 
down occurs when c > 2.8, as mentioned by Burkard 
et al£ Meanwhile, for 3D HL fails for c > 5.8, and for 
c > 1.95 for lDiiS On the other hand, at large enough R 
one always gets W > C and thus positive J, regardless 
of c. Hence, the large- R anomaly encountered in the H2 
HL solution does not occur here. We may argue that the 
large-c failure of the HL approximation is ultimately due 
to a large-^ regime [see Fig. |2Ja)]. 

The exchange calculated from Eq. J3| for the 2D sys- 
tem is given in Fig. [3] (a) for the same values of c as 
in Fig. |21 Note that for the largest c values the curves 
are interrupted when J becomes negative, and for R just 
above this point HL predicts the equally unphysical be- 
havior of J increasing with i?, 2 ° For c = 2.1, J decreases 
monotonically with R, but the "shoulder" at R ~ 4a* is 
a precursor signaling this anomaly, meaning that J > 
is not a sufficient condition for the reliability of the HL 
results. For distances far beyond these anomalous points, 
we find that the general behavior of J(R) is dominated 
by the S 2 prefactor in Eq. J2J), as shown in the inset to 
Fig. (a). 

From a formal point of view, it would be interesting to 
present a direct comparison of the R — > 00 behavior of 
J(R) calculated within HL with asymptotic expressions 
such as given in Fig. ^c) for H2 , obtained through the 
so-called median-plane method. This method was origi- 
nally introduced for two hydrogen atoms^ and has been 
extended to a variety of problems in atomic and semi- 
conductor physics. 21 It is however not directly applicable 
to the potential in Eq. © because, in the asymptotic 
limit, this model potential displays an essentially infinite 
barrier between the dots^ We present instead a quali- 
tative comparison of our results in Fig. OJa) with those 
obtained from an interpolated formula for the exchange 
coupling of two electrons under a potential related to 
ours^ namely the potential created by donor impurities 
close to a 2DEG. The connection between our results and 
those in Ref. is based on the fact that when a single 
impurity is far from the 2DEG, the electronic potential it 
creates may be approximated by a parabolic confinement, 
whose curvature defines ojq and thus c. The results plot- 
tcd in Fig.lSJb) are obtained from Eq.(44) in Ref. llSl for 
distances between the donor and the 2DEG h = 2, 4, and 
10a*, resulting in c = 2.1, 3.5 and 7 respectively. We note 
that the general trends in Figs.|3fa) and (b) are the same, 
though a detailed quantitative agreement is not found. 
This lack of agreement is due to the different nature of 
the barrier between the parabolic potentials. As noted 
above, the exchange energy within HL scales asymptoti- 
cally as S 2 , thus for the quartic model the exchange has 
an exponential decay with R 2 , while for the impurity po- 
tential in Ref. the two-electron wave- functions, and 
consequently S and J, have a slower exponential decay 
with R. 

We conclude that HL is essentially reliable, at least 
qualitatively, for the double quantum dot problem at 
large enough interdot separations. In the large over- 
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R(a*) 

FIG. 3: (Color online) (a) Exchange at long distances R, 
within HL approximation, in the two-dimensional double 
quantum dot system, for three different values of c = 2.1, 
3.5, and 7. For the larger values of c, HL breaks down at 
small R. The inset displays the curves for c = 2.1 and 3.5 
plotted versus R 2 , to show that, at large R, the decay of J 
is dominated by the 5* 2 behavior (see text), (b) Exchange as 
calculated by Ponomarev et al. (Eq. 44 in Ref. Us!) for the 
same values of c. 



lap region, HL is not applicable regardless of the sign 
of J. We have attempted a standard improvement to 
HL, which for the H2 problem is to consider the Bohr 
radius as a variational parameterii* In the present case 
this consists in taking /3 as a variational parameter: Al- 
though we obtain some lowering of the singlet and triplet 
energies by this procedure, no qualitative improvement 
is obtained. Another straightforward improvement is the 
Hund-Mulliken approach— in which the double occupancy 
states a(l)a(2) and b(l)b(2) are included in the basis set 
for the 2-electron states. This method lowers the singlet 
energy with respect to the triplet, but the failure at large 
c remains. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we have analyzed the reliability of HL 
for two types of "double-well" potentials where the 
individual-well confinement potentials arc of different na- 
tures: Atomic (Hydrogen atom) in Sec. ^ or harmonic 
(quantum dot) in Sec. II I II There are two main differ- 
ences between these problems: One refers to the strength 
of the confinement, much stronger in the first case, and 



the other to their dimensionality, 3D for atoms and 2D 
for quantum dots defined over a 2DEG. 

For the 3D H2 case, it has been long knowni^ that 
HL fails, namely, the triplet state is predicted to be- 
come the ground-state, when the atoms are very far apart 
(R ~ 50a B ). As argued above, at short distances HL is 
expected to fail as well, and it does as it overestimates 
the ground-state energy. We have also considered the 2D 
hydrogenic molecule within HL, and found qualitatively 
the same behavior and limitations as for the 3D H2 case, 
meaning that the strongly confining atomic potential is 
the dominant aspect of this problem both in 3D and in 
2D. 

Our study of the particular quartic potential pro- 
posed by Burkard et al£ for modeling coupled harmonic 
dots shows that the failure of HL at short distances is 
more dramatic here than in H2 due to the less con- 
fining potential, leading to larger overlap: A triplet 
ground-state is predicted for a wide range of parame- 
ters. We have also pointed out the fact that the quar- 
tic potential is not additive (i.e. it does not correspond 
to the superposition of the isolated harmonic potentials) 
as another possible source for this failure. Overcoming 
these limitations for the exchange calculation, in partic- 
ular when the overlap is large, requires going beyond 
the HL approximation, for example increasing the ba- 
sis set for the two- or many-electron wave-functions to 
incorporate a set of excited one-electron states. Dif- 
ferent methods have been adopted in the literature: 
Molecular orbital^^ configuration-interaction^ exact 
diagonalizationjS^SSi and unrestricted Hartree-Focki2& 
Molecular orbital calculation of double-dot model po- 
tentials have been explicitly compared to HLf2£ and in 
all cases an excellent qualitative agreement is obtained. 
According to the model potential, practically complete 
agreement may be obtained^ while deviations of up 
to 50% have also been reported^ Also, spin polarized 
DFT calculations can be done to obtain the exchange 
energy, 27 - 28 but it is difficult to estimate the reliability of 
such theories in the present scenario. Unfortunately none 
of the cited methods has the conceptual and computa- 
tional simplicity of HL, which, coupled to a careful choice 
for the model potential, should be always attempted as a 
first approximation for estimating exchange coupling in 
new systems. 
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